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Low Energy Cosmic Ray Protons from 
Nuclear Interactions of Cosmic Rays with the Interstellar Medium 

H. T. Wang 

Department of Physics and Astronomy 
University of Maryland, College Park, Maryland 20742 
The intensity of low energy (< 100 MeV) protons from nuclear 
interactions of higher energy (> 100 MeV) cosmic rays with the 
interstellar medium is calculated. The resultant intensity in the 
10 ~ 100 MeV range is larger by a factor of three to five than the 
observed proton intensity near earth. The calculated intensity from 
nuclear interactions constitutes a lower limit on the actual 
proton intensity in interstellar space. 

The knowledge of the intensity of the low energy cosmic rays 
in interstellar space is very important for the understanding of the 
origin of the cosmic rays, the heating of interstellar gas and the 
solar modulation of the cosmic rays. However, because of the expansion 
of the solar wind, cosmic ray particles suffer significant energy 
losses in the interplanetary medium, and the observed protons near 
the earth below about 100 MeV have all higher energies in interstellar 
space (Parker 1965, Goldstein y Fisk, and Ramaty 1970). The actual inter- 
stellar intensity below 100 MeV thus has to be determined by direct 
measurement on space probes at large distances from the sun. However, 
a component of this intensity can be directly calculated from presently 
available data. This component is produced by nuclear interactions 
in the interstellar medium of higher energy cosmic rays which do not 
suffer the above mentioned energy losses. 
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The purpose of the present paper is to calculate the flux of 
secondary protons in interstellar space. The nature of this calcu- 
lation is similar to previous calculations of low energy knock-on 
electrons (Abraham, Brunstein and Cline, 1966), secondary positrons 
(Ramaty and Lingenfelter , 1966), and deuterons and helium-3 nuclei 
(Ramaty and Lingenfelter, 1969). Avery crude estimate of the 
intersity of secondary protons was made by Feit and Milford (1965) 
at a time when the concept of adiabatic deceleration in the solar 
wind was not yet fully appreciated. In the present paper we use 
detailed cross sections and kinematics for proton-proton and proton- 
alpha particle interactions in order to calculate the secondary 
proton intensity below 100 MeV. 

Secondary Protons Production 

The rate of production of secondary protons by cosmic-ray inter- 
action with the interstellar material is given by: 

q s (E s ) = ^Tr^fdE j. (E) n ± a ± (E) F. (E, E g ) (1) 

t 

where is the production spectrum of secondary protons per cubic 

centimeter of interstellar space per second; E and E g are the energies 

per nucleon of the primary and secondary particles, respectively; j is 

the interstellar intensity of cosmic ray nuclei; n- L is the number of 

target nuclei per cubic centimeter of interstellar space; a is the 

interaction cross section; and F(E , E ) dE is the probability that a 

s s 

primary cosmic-ray nucleus of energy per nucleon E will produce a 

secondary proton of energy in dE around E . 

s s 
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The production is summed over all interactions i, which produce 

secondary protons. The interactions which we considered are listed 

in Table 1. We treat the secondary neutrons as protons because the 

neutrons decay to protons of essentially the same energy on a time 

scale which is negligible in comparison with the age of the cosmic 

rays. We have neglected the small contribution from the more compli- 

4 

cated interactions such as ct-He and p-CNO. It should be noted that 
a-He^ reactions are not entirely negligible but the contribution 
should not exceed 20% of the secondary proton intensity that we 
calculated . 

Both the center of mass differential cross sections for pp and 
pa elastic scattering and the total pp inelastic scattering cross 
sections were taken from Meyer (1971) . The partial cross sections for 
the different breakup modes of alpha particles in pa inelastic inter- 
actions were taken from Ramaty and Lingenfelter (1969) . The energy 
distribution F(E,E g ) for pp and pa elastic scattering are directly 
obtained from the differential cross sections. The distribution 
F(E,E s ) for pp inelastic scattering is taken from the experimental 
data compiled by D. Lai (private communication). The method for the 
evaluation of F(E,E g ) for pa inelastic interactions is similar to that 
of Ramaty and Lingenfelter (1969) for deuterium and helium-3 production. 

A detailed discussion of cross sections and kinematics used in the 
present calculation will be given in a forthcoming thesis (Wang, University 
of Maryland, dissertation) . 

For the primary proton and alpha particle spectra in interstellar 
space, we use power laws in total energy with low energy cutoffs at 
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100 MeV per nucleon. Such spectra were found to be consistent 
with the modulation of both protons and electrons in the inter- 
planetary medium and the non-thermal radio background produced by cosmic 
electrons (Goldstein, Ramaty, and Fisk 1970) . 

Using equation (1) , we can now calculate the production rate of 
the secondary protons. The production of each breakup mode for pa 
and ap inelastic interactions are shown in Figure 1 and Figure 2, 
respectively. We show in Figure 3 the production spectra of pp elastic 
and inelastic scattering, pa and ap elastic scattering, and the sum 
of pa and ap inelastic interactions. The low energy cutoff in the pa 
elastic scattering production spectrum results from the low energy cut- 
off at 100 MeV per nucleon and the fact that it is impossible for a 
light particle to give all its energy to a heavier target particle in 
one interaction. The complicated structure of the spectra in Figure 1 
and Figure 2 are the results of the following reasons. First we have 
added the spectra with pion production to that without pion production. 
Second we have evaluated the kinematic quantities of multinucleon particles 
(d, He J etc) first and these particles have a minimum secondary energy 
about 10 MeV per nucleon in ap inelastic interactions (Ramaty and 
Lingenfelter , 1969) and hence a same cutoff in proton production is 
obtained. Third, as in ap elastic scattering, there is a cutoff in the 
primary intensity. We note that ap plus pa inelastic interactions con- 
tribute approximately 30 percent of the total production. Because several 
assumptions have been made about the kinematics of these interactions, 
an uncertainty of thirty percent is introduced. 
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Secondary Proton Intensity in Interstellar Space 
From the production spectrum in interstellar space computed 
above, we may determine the equilibrium density of secondary protons 
in the galaxy. We use the following equation: 


N„ , 3 dE 

7 <». 


( 2 ) 


where N g is the secondary equilibrium density; q g is the secondary 

production rate calculated above; t is the effective leakage life 
dE 

time, is the rate of energy loss due to ionization. Changing 
parameter from time t to matter traversal x = pvt, where J> is the density 


„-3 


of interstellar material in gem , v is the velocity of the secondary 


particle, we find that the secondary intensity j is given by: 


j s (E) = 77 L dE' q (E ’) exp (- i- f E dE” ) (3) 

4ttP iEJ E s X 0 J E dE^ 

' dx dx 

wHere m p is the proton mass and x Q = j?vx. The numerical value for x D is 
about 5g cnT^ (Ramaty and Lingenfelter , 1971). Below 500 MeV, the 
ionization loss rate ip- in MeV g"^cm^ can be approximated by 102 *® 

(Barkas and Berger, 1964). The resultant proton intensity is shown in 
Figure 4.. Also shown in Figure 4 are the primary proton interstellar 
intensity that produces the secondary protons by the nuclear interactions. 
The data points represent the observed intensity near earth (Fan, et al, 
1966; Balasubrahmanyan, et al., 1966 (a, b) ; Wadding ton and Freier, 1966; 
Ormes and Webber, 1966; McDonald, 1968). As can be seen the calculated 
secondary intensity is larger by a factor of three to five than the observed 
intensity and could be measurable by detectors on space probes at large 


distances from earth. 
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In order to get some idea about the magnitude of this intensity 
in comparison with the low energy proton intensity that may exist 
in interstellar space, we calculate the rate of ionization £ of the 
interstellar medium by the s econdary intensity shown in Figure 4. 

We find that £ = 3 x 10 (sec Hatom) . This value is several order 
of magnitude lower than values of |(10"^~10 - ^) obtained from Ha and 
Hg observations (Reynold, Roesler, and Sherb, 1972). However, because 
the ionization of interstellar medium is not necessarily due to low 
energy protons, the observed values of £ only set upper limits on 
the low energy proton intensity in interstellar space. The present 
calculation constitutes a lower limit on the intensity of these particles. 

Acknowledgement: I wish to thank Dr. Reuven Ramaty for valuable 
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Table 1 


Secondary Proton Production 


1. 

P 

4- H 1 

~y 

p + p (elastic) 

2. 

P 

+ H 1 

-> 

p + p + (7 t) (inelastic) 

3. 

P 

+ He 4 

-> 

p + a (elastic) 

4. 

a 

+ H 1 

->■ 

p + a (elastic) 

5. 

P 

+ He 4 

-* 

t + 2p + (tt) 

6. 

a 

+ H 1 


t + 2p + (tt) 

7. 

P 

+ He 4 

-► 

He ^ + p+ n+ (it) 

8. 

a 

+ H 1 

-> 

He^ + p + n 4- (it) 

9. 

P 

+ He 4 


d + 2p + n + (tt) 

10. 

a 

+ H 1 

-> 

d 4- 2p 4- n 4- (tt) 

11. 

P 

+ He 4 

-f 

2d + p + (tt) 

12. 

a 

+ H 1 


2d + p + ( 7r ) 

13. 

P 

+ He 4 

-* 

3p + 2n + (tt) 

14. 

a 

+ Hi 

-* 

3p + 2n + (tt) 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Proton production spectra for different breakup modes 
for po< inelastic reactions shown in Table 1. 

Proton production spectra for different breakup modes 
f or o^p inelastic reactions shown in Table 1. 

Proton production spectra for pp elastic and inelastic 
scatterings, 0(p and pe^ elastic scatterings, and pc< 
plus o<p inelastic reactions. 

Proton intensities. 
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